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Abstract 

Several extended supersymmetric models, motivated by either grand unification, or by 
neutrino mass generation, predict light doubly charged higgsinos. We study the produc- 
tion of a single doubly charged higgsino and its decay channels at the International Linear 
Collider (ILC) operating in the e~e~ mode. We analyze the production cross section for 
e~e- -> A7 R X\ as a function of different kinematic variables, followed by the decay, through 
several channels, of the doubly charged higgsino into a final state of two leptons plus missing 
energy. We include the standard model background and discuss how kinematic cuts could 
be used effectively to limit this background. Single production of these exotics could provide 
a spectacular signal for a new underlying symmetry and for physics beyond the minimal 
supersymmetric standard model. 
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1 Introduction 



It is widely expected that the standard model (SM), though successful at predicting almost 
all experimental data in low energy physics, is not the complete theory of fundamental inter- 
actions. In addition to the theoretical inconsistencies and incompleteness of the theory, there 
have been recently experimental incentives to study models beyond the SM. Observations 
and measurements of solar and atmospheric neutrino oscillations [1], as well as indications 
of hot and cold dark matter [2] are not predicted, and cannot be explained, respectively, by 
the SM. 

Supersymmetry, in the form of the minimal supersymmetric standard model (MSSM) 
is the most popular scenario of physics beyond the SM. It provides a satisfactory (at least 
qualitatively) explanation for dark matter, but suffers from the same problems as the SM 
when it comes to explaining neutrino masses. One must either invoke R-parity violation [3] 
and abandon the dark matter candidate, or add right-handed neutrinos and introduce the 
see-saw mechanism [4] . Supersymmetric grand unified theories (SUSY GUTs) which contain 
left-right supersymmetry resolve both problems naturally [5,6]. An added bonus of SUSY 
GUTs is that electromagnetic, weak and strong interactions unify at the same energy scale. 
If such SUSY GUTs containing left-right supersymmetry are present in nature, one must 
devise methods to search for signals inherent in them, and absent in other models. One such 
signal would be the discovery of exotic particles unique to such models. 

Supersymmetric left-right theories (LRSUSY), based on the product group SU(2) L x 
SU(2) R x U(1)b-l, are attractive for many reasons [5,7]. They disallow explicit R-parity 
breaking in the Lagrangian; they provide a natural mechanism for generating neutrino 
masses; and they provide a solution to the strong and weak CP problem in MSSM [8]. 
Neutrino masses are induced by the see-saw mechanism through the introduction of Higgs 
triplet fields which transform as the adjoint of the SU (2) R group and have quantum numbers 
B — L = ±2 (where B is baryon, and L is lepton number). While the Higgs triplet bosons 
are present in the non-supersymmetric version of the theory, their fermionic partners, the 
higgsinos, are specific to the supersymmetric version. It has been shown that, if the scale for 
left-right symmetry breaking is chosen so that the light neutrinos have the experimentally 
expected masses, these higgsinos can be light, with masses in the range of (9(100) GeV [9]. 
Such particles could be produced in abundance at future colliders and thus give definite signs 
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of left-right symmetry at future colliders. 

While the production of doubly charged Higgs bosons in the framework of the left-right 
model has been investigated at linear accelerators [10] and at LHC [11], the corresponding 
higgsinos have received less attention. The exceptions are the references in [12], where some 
of the properties of doubly charged higgsinos have been highlighted, and in [13], where the 
pair production of higgsinos in e + e~ was analyzed. We note that doubly charged higgsinos 
also appear in the so-called 3-3-1 models (models based on the SU(3) C x SU(3)l x U(1)n 
symmetry) [14]. 

In this present work, we concentrate on a definite signal for doubly charged higgsinos: 
the production of a single one at an e~e~ collider, followed by the decay (through several 
channels) to e~e~ + E miss . In order to obtain definite predictions for the signal, we perform 
our analysis in the context of the LRSUSY, though we expect the results for the 3-3-1 model 
to be similar. We concentrate first on the details and characteristics of the production 
cross section. We then discuss the possible decay modes of the doubly charged higgsinos 
(either two-body or three-body, depending on the spectrum characteristics). We complete 
our analysis with the discussion of the SM background, and indicate how cuts could be 
employed most efficiently to reduce these backgrounds. 

Our paper is organized as follows. In Section [2] we describe the LRSUSY model, with 
particular emphasis on the sectors of interest, as well as summarize the restrictions on the 
relevant Yukawa couplings. In Section[3]we present the analysis for the production and decay, 
separately for the left and right-handed doubly charged higgsino and of the SM backgrounds 
and cuts needed to observe the signal. We reach the conclusions in Section HI The Appendix 
contains the mixing in the scalar and gaugino/higgsino sectors which enter our calculation. 



2 Description of the LRSUSY Model 



The minimal supersymmetric left-right model is based on the gauge group SU (3)c x SU{2) l x 
SU(2)r x U{V)b-l- The matter fields of this model consist of three families of quark and 
lepton chiral superfields which transform under the gauge group as: 
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where the numbers in the brackets denote the quantum numbers under SU(3)c x SU{2) L x 
SU(2) R x U{1)b-l- The Higgs sector consists of the bidoublet and triplet Higgs superfields: 
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The bi-doublet Higgs superfields $i,$2 break the SU(2) L x C/(l)y symmetry and generate 
a Cabibbo-Kobayashi-Maskawa mixing matrix. Supplementary Higgs representations are 
needed to break left-right symmetry spontaneously: triplet Higgs Ax,, A r bosons are chosen 
to also support the seesaw mechanism. Since the theory is super symmetric, additional triplet 
superfields Sl, 8r are needed to cancel triangle gauge anomalies in the fermionic sector. The 
most general superpotential involving these superfields is: 

W = Y^Q T ^r 2 Q c + Y^ ) L T ^ i ir 2 L c + i(h n L T T 2 5 L L + h n L cT r 2 A R L^ 

+/i 3 [Tr(A L 5 L + A R 8 R )} + mjTrfaQjiTzQj) + W NR (3) 

where Wnr denotes (possible) non-renormalizable terms arising from higher scale physics or 
Planck scale effects [9]. The presence of these terms insures that, when the SUSY breaking 
scale is above Mw R , the ground state is R-parity conserving. In addition, the potential also 
includes F-terms, D-terms as well as soft supersymmetry breaking terms: 

'A^Y^Q T ^tT 2 Q c + A i I Yf~L T ^ l iT 2 L c + iA LR h u (L T r 2 5 L L + L cT r 2 A R L c ) 
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The symmetry is broken spontaneously to U(l) em when the neutral Higgs fields acquire 
non-zero vacuum expectation values (VEV's): 
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The Lagrangians in ([3]) and (jl]) give rise to the interactions of the doubly charged higgsinos 
A~[~ R . As these are determined by the magnitude of the triplet Yukawa couplings, we review 
first the restrictions on these. 



2.1 Experimental limits on triplet masses and Yukawa couplings 

Indirect experimental limits for the triplet Yukawa couplings come from lepton number 
violating processes mediated by doubly charged Higgs bosons. These processes constrain 
the first two generations, while the third generation couplings remain constrained only by 
the requirement of perturbativity. The constraints on Yukawa couplings as a function of the 
doubly charged Higgs mass can be found in [15,16] and are given as follows 

h efl h ee < 3.2 x 1CT 11 GeV~ 2 • M A — from/x -> eee, 
h^h^ < 2 x 1(T 10 GeV~ 2 • M A __ from/x -> ej, 
h\ e < 9.7 x 1(T 6 GeV~ 2 • M\— from Bhabha scattering, 
h 2 ^ < 2.5 x 1(T 5 GeV" 2 • M A __ from (g - 2) M , 

/i ee V < 2.0 x 1(T 7 GeV~ 2 ■ Ml— from muonium — antimuonium transition. 

Bhabha scattering indirect limits were studied also at LEP, and it was found that if 
h ee > 0.7, the doubly charged Higgs mass should be in the TeV region [17, 18]. Also at 
LEP, direct searches of doubly charged Higgs were performed. From the pair production 
of doubly charged Higgses a lower bound around 100 GeV was established, if hij > 10~ 7 , 
i,j = e,fi,r [18,19]. From the single production of doubly charged Higgs boson in e + e~ 
collisions, couplings h ee < 0.071 are allowed, if Ma — < 160 GeV, assuming 100% branching 
fraction to leptons [17]. At HERA the single production of doubly charged Higgs boson was 
studied in e + p collisions [20] and a lower limit of 141 GeV was found on the doubly charged 
Higgs mass if the coupling h efJ , = 0.3, while the mass limit for h eT = 0.3 was 112 GeV. For 
heavier doubly charged Higgses, the corresponding Yukawa couplings are less restricted, and 
for Ma — > 150 GeV triplet Yukawa couplings are not constrained by HERA. Note that the 
direct search limits are the only ones which include the third generation couplings. 
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For simplicity we assume that the flavor-violating off-diagonal entries in the coupling 
matrix are zero and only the flavor diagonal elements have finite entries. Assuming that 
the triplet Yukawa couplings are degenerate the most stringent constraints for the AL = 2 
couplings come from results for the muonium-antimuonium transitions which depend on the 
mass scale for the triplet Higgs fields viz. h u = f u < 0.44M A ±± TeV~\ at 90% C.L. [16]. 

We now proceed to examine the doubly charged higgsino sector of the theory. We first 
review the known information on their masses and couplings, then proceed to discuss their 
production and decay channels. For completeness, we present the information on the mixing 
in the chargino, neutralino and charged scalar sectors of our theory, insofar as relevant to 
our discussion of doubly charged higgsinos, in the Appendix. 



3 Production and Decay of Doubly Charged Higgsinos 

We analyze first the single production of the doubly charged higgsino at the next generation 
linear collider running in the e~e~ mode. The linear collider operating in this mode will 
provide an ideal environment for single production of such doubly charged particles and 



e A 




Figure 1: Feynman graphs for the two-body production of A ^{ at e e collider. 



signals for such processes can be remarkably free from backgrounds. We focus on a typical 
production mode for the doubly charged higgsinos A in association with the lightest 
neutralino and calculate the rates for the process: 

e~e- — A~x° (5) 
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and study the signal resulting from the decays of the A . This production is mediated 
by the selectron exchange^ as shown in Fig [1] and would be accompanied by large missing 
energy as the lightest neutralino in the final state is stable and escapes undetected. The 
e~e~ production mode would allow us to probe a large range of masses of the doubly charged 
higgsino, compared to the case where one produces them in pairs at linear colliders operating 
in the e + e~ mode. The upper value of the mass of the produced A will be constrained only 



Sample point A : tan/3 = 20 

M B _ L = 25 GeV M L = M R = 250 GeV 

[L X = 1000 GeV fj, 3 = 300 GeV 

v Al = 1.5 x 10~ 8 GeV v Sl = 1.0 x 10~ 8 GeV 

VAr = 3000 GeV v Sr = 1000 GeV 


M~o = 91.8 GeV M~ = 180.6 GeV 

Al A 2 

M~o = 250.0 GeV 

M~± = 249.0 GeV M~± = 300.0 GeV 

Al A2 

M~± = 911.7 GeV 

A3 


Sample point B : tan/3 = 30 

Mb-l = 100 GeV M L = M R = 500 GeV 

/il = 500 GeV /i 3 = 500 GeV 

v Al = 1.5 x 10~ 8 GeV v 5l = 1.0 x lO" 8 GeV 

v Ar = 2500 GeV v 5r = 1500 GeV 


M~o = 217.3 GeV M~ = 441.7 GeV 

Ai A2 

M~o = 450.0 GeV 

A3 

M~± = 447.8 GeV M~± = 500.0 GeV 

Al A2 

M~± = 500 GeV 

A3 



Table 1: Sample points of the particle spectrum. The left column contains the input values for 
the parameters of the model, while the right column lists the masses of the lightest the neutralinos 
and charginos corresponding to the choice of input parameters. 



by the mass of the lightest neutralino and the kinematic threshold of the machine energy 
of the collider. In the latter part of this section we will focus on the possible signatures for 
different final states. 

For the analysis, we choose two representative points in the parameter space of the model 
as presented in Table [T] (we have consistently chosen negligibly small values for the VEV's 
va l and v$ L so that the constraints on the p-parameter is not disturbed). This will lead, as 
our discussions will show, to clear signals for doubly charged higgsino production. Moreover, 



1 We neglect the contribution coming from the exchange of doubly charged Higgs (A) in the s-channel 
which we find to be negligibly small in comparison. This is because of the huge suppression coming from 
the mass of A in the propagator, which we have assumed to be large. 
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Figure 2: The production cross section as a function of the triplet higgsino mass for y/s = 1 TeV. 
We highlight the dependence of the production cross section (for A^~J on the beam polarization 
(-1,-1) as compared with the unpolarized cross sections for two choices of selectron mass (250 GeV 
and 500 GeV). The plot is produced for the sample point A given in TabieQJ 




a reasonable choice for the AL = 2 couplings leads to like-sign lepton signals which can 
have very suppressed or no SM background. For simplicity, we have however focused only 
on the diagonal AL = 2 couplings. As one could produce both the left or right-chiral states 
of A , we divide this section into two parts and discuss the production of the two chiral 
states separately. 

3.1 The left-handed higgsino A x 
3.1.1 Production 

In Fig.[2]we plot the production cross sections of A^~ with both unpolarized and polarized e~ 
beams. The center-of-mass energy is taken to be y/s = 1 TeVo We choose a coupling strength 
of f ee = 0.1, in agreement with bounds on AL = 2 couplings from different experimental 
data as listed in section 12.11 With masses for the triplet Higgs fields in the TeV-range one 



2 While the initial center-of-mass energy at the ILC [21] is expected to be y/s = 500 GeV, it is likely to 
be increased soon to 1 TeV, for which more promising signals can be observed. 
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can easily obtain larger values for the couplings. In contrast to existing bounds on the mass 
of the triplet Higgs (quite severe for larger values of the AL = 2 couplings), the mass of A^~ 
is not strongly constrained (<; 100 GeV). So a light with a large AL = 2 coupling is 
not ruled out and remains consistent with experimental data. In the sample points given in 
Tabled] the lightest neutralino is the lightest supersymmetric particle (LSP). We have used 
sample point A in calculating the rates in Fig. [2j 

The production cross sections are shown for two different choices of selectron mass viz. 
250 GeV and 500 GeV. As expected, the larger value of the selectron mass in the propagator 
suppresses the production cross section. The production is however enhanced by a factor 
of four for the left-chiral A^~, if one uses 100% left-polarization for both beams, as shown 
in Fig [2j Other combinations make the production cross section vanish for the left-chiral 
state because of the (1 — 75) coupling (see Section [53]) . The figure shows that even with a 
conservative value for the AL = 2 coupling, one expects a production cross section between 
10-14 fb for polarized beams and a 250 GeV-mass selectron. With an integrated luminosity 
of 500 fb~ l one can have a large production rate of the doubly charged higgsino Aj^ . 

It is worth noting that the choice of the soft parameters have a role to play in the pro- 
duction mechanism, as they form the basic entries in the matrix diagonalizing the neutralino 
mass matrix (see Section [5721 . The e — e — x° coupling, which can be written down as 

rllxl - ^(9LN kl + g v N k3 )P L 

and the components of the lightest neutralino mass eigenstate will thus be modified for 
different choices of the soft parameters in the theory. A quick look at Table [1] shows that 
the mass splitting between the lightest neutralino and the next lightest neutralino strongly 
depends on the ratio where a larger ratio results in smaller mass splitting. This can 

VS R 

lead to interesting signatures where the next lightest neutralino can become the next-lightest 
SUSY particle (NLSP). Such interesting spectrum analysis is however left for future work 
and not considered here anymore. We show the dependence of the relevant soft parameters 
on the lightest neutralino mass and the corresponding production cross sections in Fig. [3] 
and Fig. H] The mass of the doubly-charged higgsino is Mr = U3, with the exchanged 

L 

selectron mass fixed at mg L = 500 GeV while the polarization choice for the colliding electron 
beams is (-1,-1) which maximizes the production cross section for the left-handed higgsino. 
In Fig. [3](a), we show the lightest neutralino mass as a function of the U{1)b-l gaugino 
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Figure 3: (a) The lightest neutralino mass as a function of Mb-l- (b) The production cross section 
as a function of Mb-l, for three different values of fi 3 . The other parameters are matched with 
the remaining inputs to sample point A given in TaWed Here yfs = 1 TeV, rrie L = 500 GeV and 
Af A -- = /13. 

mass parameter Mb-l for three different values of pL 3 . We find that, increasing value of 
Mb-l (the other model parameters fixed to values for sample point A given in Tabled]), the 
lightest neutralino state becomes heavier. In fact the the lightest neutralino mass shows no 
dependence on /j, 3 for low values of Mb-l, but can be seen to admit a dominant admixture of 
the right higgsino state beyond Mb-l ~140 GeV for the lower choice of fi 3 = 300 GeV. This 
fact is highlighted in the plot for the cross section, which starts falling rapidly for values of 
Mb-l > 140 GeV as shown in Fig.[3](b), for /i 3 = 300 GeV. For large Mb-l, the lightest state 
contains a dominant admixture of the right higgsino states and thus its coupling strength 
is decreased. The cross section for the other choices of fi 3 = 500 GeV and 700 GeV do not 
show the sharp fall as the lightest neutralino remains dominantly U(1)b-l~ gaugino-like for 
the choice of M B -l in the plots. 

In Fig. Hl(a) we plot the lightest neutralino mass as a function of Ml = Mb = M 2 for 
the three different values of li 3 . The mass of \i again shows little dependence on /i 3 for low 
values of M 2 . The contribution of the right higgsino state in the lightest neutralino shows 
up early for the lower values of /13 as seen in Fig. |U(a). As shown in Fig. IU(b), in contrast 
to Fig. [3](b), the cross section varies over a larger value when increasing the value of the 
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Figure 4: (a) The lightest neutralino mass as a function of Ml = Mr = Mi- (b) The production 
cross section as a function of Mi , for three different values of . The other parameters are matched 
with the remaining inputs to sample point A given in TaWeHJ Here = 1 TeV, ra^ L = 500 GeV 
and — = ^3. 

neutralino mass (increasing value of the soft parameter). The cross section is again seen to 
show a sharp fall at different values of Mi for different choices of /X3. This is again due to the 
increase of the right- higgsino component in the lightest neutralino and the corresponding 
change in the values of the e — e L — X\ coupling. The curves for /x 3 = 700 GeV end abruptly, 
because the cross section becomes zero as the lightest neutralino mass becomes greater than 
300 GeV and the process is kinematically disallowed at the \fs = 1 TeV machine. 
We now proceed to perform the signal analysis for the production in detail. 

3.1.2 Decays 

We focus on the event rates for specific final states arising through the decay of the left- 
chiral doubly charged higgsino. This would require the knowledge of all possible decay 
modes accessible to the A . In general, the following 2-body decays of the doubly charged 
higgsinos are allowed: 

• A — >£- £~, 

• A — A" xl 
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• A — > xT A". 

However, we assume the triplet Higgses to be heavier (which allows a less constrained 
AL = 2 coupling), and degenerate in mass, which renders them kinematically inaccessible 
for decay modes of the relatively lighter doubly charged higgsinos. Hence, we concentrate 

on the only favored channel for decay which is A — > t~ £~, provided mj < . For 

relatively light higgsinos, one can in principle have mj > and in such a scenario the 

only allowed decay mode would be the 3-body decays for the doubly charged higgsinos, which 
would dominantly be through off-shell sleptons: A — > i*~ £~ —>■ £~t~%\. We focus our 
analysis on the final states: two like signed leptons in association with large missing energy 
coming from the LSP (xT)- 

3.1.3 Analysis of the final states 

In the rest of the analysis of final states, we will consider the 2-body and 3-body decay modes 
of the doubly charged higgsino (A^~), and look at the resulting signal events against the 
most dominant SM background. We present our results for both the sample points given 
in Table Q] with point A used for a machine with a/s = 500 GeV and point B considered 
for a/s = 1 TeV option. In the case of sample point A, the mass of the doubly charged 
higgsino (A^ - ) is taken to be 300 GeV and the LSP mass is 91.8 GeV, while the other 
neutralino states are heavy and not accessible to the decay of the doubly charged higgsino 
or sleptons (directly). This also makes both the 2-body A — ► £r £~ and 3-body decay of 
A — > £~£~Xi the only allowed channels of decay. However, we point out that if the other 
decay modes > 2 and xt u i) were accessible as the decay channels of sleptons, then 

the only visible effect on our signal would be to reduce the branching fractions and hence 
decrease the signal events for our final states. But they will not change the characteristic 
kinematic features of the final states in consideration. In our analysis we focus on the 
following leptonic final states: 

(i) e-e-fi (ii) yTyr$ 

It is worth pointing out that the analysis for a final state t~t~ f, will be exactly similar to 
pT [flfj, provided that the AL couplings are identical and that we assume similar detection 
efficiencies for the r's in the final state. Thus we do not discuss the final states with r's 
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in our analysis any further. Assuming, as before, beam polarization of (-1,-1) throughout 
the analysis for the A~[~ production, the major SM background that contributes to the final 
states in (i) is the scattering process: 

e" + e~ — > e" + e~ + V{Vi 

which, although a continuum background, could prima jade be very large. Thus, the event 





^ = 500 GeV 


Cuts Used 


SM 


signal- 1 


signal-2 


E e > 5 GeV |r/ e | < 3.0 
|;>10GeV AR ee > 0.2 


(-,-) 537.7 fb 
(+,+) 15.1 fb 


(-,-) 123.8 fb 
(+,+) 486.8 fb 


(-,-) 19.9 fb 
(+,+) 78.1 fb 


E e > 5 GeV \r] e \ < 1.5 
|;>100GeV Ai? ee > 0.2 


(-,-) 218.0 fb 
(+,+) 3.8 fb 


(-,-) 102.7 fb 
(+,+) 403.98 fb 


(-,-) 16.5 fb 
(+,+) 65.1 fb 


E e > 5 GeV |r/ e | < 3.0 
^>^/2GeV Ai? ee >0.2 


(-,-) 280.1 fb 
(+,+) 3.3 fb 


(-,-) 123.8 fb 
(+,+) 468.8 fb 


(-,-) 19.9 fb 
(+,+) 78.1 fb 


E e > 5 GeV |r/ e | < 1.5 
^>Vs/2GeV Ai? ee >0.2 


(-,-) 103.8 fb 
(+,+) 0.103 fb 


(-,-) 102.7 fb 
(+,+) 403.98 fb 


(-,-) 16.5 fb 
(+,+) 65.1 fb 




^ = 1 TeV 


E e > 5 GeV |r/ e | < 3.0 
|;>10GeV A^ > 0.2 


(-,-) 1.13 pb 
(+,+) 12.6 fb 


(-,-) 40.5 fb 
(+,+) 156.4 fb 


(-,-) 14.0 fb 
(+,+) 53.9 fb 


£ e > 5 GeV \ri e \ < 1.5 
#> 100 GeV AR ee >0.2 


(-,-) 238.9 fb 
(+,+) 3.1 fb 


(-,-) 33.4 fb 
(+,+) 129.2 fb 


(-,-) 11.7 fb 
(+,+) 45.0 fb 


E e >5 GeV |r/ e | < 3.0 
E>^/2GeV AR ee >0.2 


(-,-) 605.9 fb 
(+,+) 0.4 fb 


(-,-) 40.5 fb 
(+,+) 156.4 fb 


(-,-) 14.0 fb 
(+,+) 53.9 fb 


E e >5 GeV \rj e \ < 1.5 
|> V«/2 GeV Ai? ee > 0.2 


(-,-) 106.0 fb 
(+,+) 0.007 fb 


(-,-) 33.4 fb 
(+,+) 129.2 fb 


(-,-) 11.7 fb 
(+,+) 45.0 fb 



Table 2: Signal and SM cross sections for the e~e~$ final states with different choice of kinematic 
cuts for both signal and background at the e~e~ collider with center-of-mass energies y/s = 500 
GeV and yfs = 1 TeV. We also show the beam polarizations in parentheses. The AR is defined 
as (AR) 2 = (Acf)) 2 + (At/) 2 with An and Acf) respectively denoting the separation in rapidity and 
azimuthal angle for the pair of particles under consideration. 

selection criteria are largely aimed at suppressing this continuum background. The SM 
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background for signal (ii) comes from a six-body final state 

e" + e~ — > pT + \i~ + z/ M z/ M + z/ e i/ e 

and is quite small. 

To highlight the effect of kinematic cuts in suppressing the background, we present some 
results in Table [2j The SM background has been calculated using the event generation 
package of Madgraph and Madevent [22] , with slight modifications to extract the polarized 
cross sections. We have chosen the |AL| = 2 coupling strength as f ee = 0.3 throughout the 
analysis and the branching ratios for the 2-body and 3-body decays as: 

BR(A L - -> £- L £-) =-, i = e,fi, r, m~ la < M K -- 
BR(£; L ^ £;xi) = l, i = e,/i,T 
BR(A L - - tre-jg) = i i = e,n, r, m~ la > M A -- 

where only the 3-body decay is allowed when > To define our notations, we 

write signal- 1 to correspond to the 2-body decay of A^~ and write signal- 2 to correspond 
to the 3-body decay of A^~. We use sample point A for the machine with y/s = 500 GeV 
as its center-of-mass energy and sample point B for the y/s = 1 TeV option. We choose 

— = 300 GeV and m^. L = 150 GeV for signal- 1, while for signal-2, = 400 GeV to 
study the signal at the y/s = 500 GeV machine. The corresponding choice for the analysis 
at the y/s = 1 TeV option is — = 500 GeV, m|. £ = 250 GeV for signal-1, while for 
signal-2, m- t = 550 GeV. 

The background for the signal (ii) fi~fi~f! is a six-body final state and would be very 
small. The cuts, which gave the SM background a(e~e~ $)= 537.7 fb/1.13 pb, give cross 
sections: cr(fi~fi~f) ~ O(10 _1 )/O(l) fb for, respectively, the = 500/1000 GeV machine. 
However the effect of cuts are not severe for the n~fi~$ final states as compared to e~e~]fi 
final states, where the cross section is reduced by more than ~ 80 — 90% of the original one, 
while the cross section in this case reduces by about 30% for the more stringent cuts listed 
in Table [2j Hence we have chosen to use the same set of cuts for both final states. 

We present our results for the set of cuts which reduces the background while keeping the 
signal relatively large. Our choice of cuts for the signal-background analysis is as follows: 

• We demand that the electrons respect a minimum rapidity cut of \q e \ < 1.5 which is 
the most effective cut to reduce the large continuum background. The signal is peaked 
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at r] e = 0, while the SM background is peaked for \rj e \ > 1.5, which is predominantly 
due to the strong t-channel photon contribution in the e~e~fi final states. 

• The final state electrons must carry a minimum energy, E e > 5 GeV. 

• We also demand a large cut on the missing energy because of the massive LSP in the 
final state for the signal, ]fi> \fs/2 GeV. 

• To ensure proper resolution between the final state electrons we demand that they are 
well separated in space and satisfy AR ee > 0.2. 

Using the above cuts, we can see from Table [2]that the large SM background which previously 
overwhelmed the signal for the most conservative cuts, is effectively suppressed without 
loosing out much on the signal events. 

We present the different kinematic distributions in figures [51 El [TJ [8] and [9] for the final 
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Figure 5: Binwise distribution of missing energy for both signal-1 and signal-2 for the final states 
e~e~ (y/s = 500 GeV). The broken (— ■ ■—) brown lines correspond to the signal through 
production while the broken (—■■■—) green lines correspond to the signal events from e^-pair 
production. Also shown in dashed and dotted dark lines are the 5a and 3a fluctuations in the SM 
background. Each binsize is 10 GeV. (a) £ = 100/6" 1 , (b) £ = 500/6" 1 . 



states e e J5for the 500 GeV machine. We must point out that the SM background for the 
other final state \T ' n,~fta,t the 500 GeV machine turns out to be completely insignificant and 
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the signal would be background free. In fact the signal would also act as an obvious hint for 
scenarios with doubly charged higgsinos and provide strong distinguishing features from the 
typical MSSM signal coming from selectron pair production (e~e~ — > e^e^), which produces 
an identical final state e~e~^but gives no contribution to the fi~fi~fi final states. Having 
said this, we must note that for the e~e~fi final states, we have additional events coming 
from e~e~ — > e^e^ in this model too. Looking at our mass spectrum, this production mode 
will only contribute to signal-1. We estimate this to be 29.72 fb at the 500 GeV machine 
for the polarization choice (-1,-1) and the same set of cuts. We include this contribution in 
signal-1 for all the figures ([5H9]), and the integrated luminosity taken in (a) and (b) is 100 
and 500 fb^ 1 respectively. We have also shown the independent contributions for signal-1 
coming from and e^-pair production in Figs. [5], [6] and [9] to highlight the difference in the 
two cases. As mentioned earlier, the doubly charged higgsino mass is taken to be 300 GeV, 
while the selectron masses are 150 GeV and 400 GeV for signal-1 and signal-2 respectively. 
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Figure 6: Binwise distribution of the invariant mass of the visible particles in the final state e~e~f] 
for both signal-1 and signal-2 (yfs = 500 GeV). The broken (— ■ ■—) brown lines correspond to 
the signal through production while the broken (—•"—) green lines correspond to the signal 
events from e~[-pair production. The background follows the notation of Fig\E[ Each binsize is 10 
GeV. (a) C = 100 fb' 1 , (b) £ = 500 fb' 1 . 



In Fig. [5] we plot the binwise distribution of missing energy for both signal-1 and signal-2, 
represented by colored histograms, and compare it with the expected statistical (Gaussian) 
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fluctuations at the 3 and 5 standard deviations of the SM background. As expected, the 
signal would be associated with a large missing energy due to the presence of LSP in the final 
state, and this also allows to set a strong cut on the missing energy. In Fig. [5(a) we show 
the distribution for an integrated luminosity of 100 fb -1 and find that signal- 1 stands out 
quite distinctly against a 5a fluctuation in the SM background. The fact that the signal is 
actually superimposed over the tail end of the much wider missing energy distribution of the 
SM, carried by the neutrinos, makes signal-1 stand out against such large fluctuations in the 
SM background. We show the same distribution in Fig. [5(b) with an integrated luminosity 
of 500 fb" 1 , where both signal-1 and signal-2 stand out against a 5cr fluctuations in the 
background. 

In Fig. [6] we show the invariant mass distribution for the pair of visible particles in the final 
state. As both the visible leptons come from the (cascade) decay of A^~, the distribution for 
the signal would depend strongly on the masses of and x?' s i n the final state, and the 
distribution exhibits a sharp kinematic edge which highlights this dependence. The location 
of this kinematic edge can be well approximated by the formula 

5oooi 1 1 1 1 11 11 1 1 11 1 1 11 11 1 11 11 1 1 11 1 ioooo 
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Figure 7: Binwise distribution of the rapidity of the e in the final state e e ft for both signal-1 
and signal-2 (y/s = 1 TeV). The background follows the notation of Fig\^ Each binsize is 0.2 
radians, (a) £ = 100/Jr 1 , (b) C = 500/6" 1 . 
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where we have assumed that the final neutralino coming from the cascade decay of A~[~ 
carries a negligible 3-momentum, as its mass is nearly a third of the mass of A^~. This is in 
sharp contrast to the SM distribution, where the background has a more uniform invariant 
mass distribution. In fact, the distribution arising from the independent contribution of 
e^-pair production shown in broken (—•••—) green line has a uniform distribution like the 
SM background, but is overwhelmed by the events coming from the A^~ production. The 
luminosity choices are as for the plots in Fig. [6](a) and (b). 
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Figure 8: Binwise distribution of the AR between the visible particles in the hnal state e e ft for 
both signal- 1 and signal-2 (yfs = 500 GeV). The background follows the notation of Fig Each 
binsize is 0.2. (a) C = 100/6 -1 , (b) C = 500/6 -1 . 



In Fig. [7] and Fig. [8], we give the rapidity distribution for the final state leptons and 
the AR distribution between the pair of visible particles in the final state respectively. The 
rapidity distribution clearly suggests that the leptons in the signal distribution emerge back 
to back, since the distribution peaks at rj = 0. This is because the leptons are decay products 
of a massive A which has low boost, and hence it is more likely to give off back to back 
lepton-slepton in signal- 1, while the other lepton coming from the heavy slepton is most 
likely to carry the maximum boost of the slepton as it is produced against a massive particle 
(x?) with less boost. The distribution for signal-2 shows a similar behavior because the 
2-lepton system recoils against a massive LSP which would have a minimum boost. For 
the background, the distribution is mostly peaked at large values of \rj\ due to the strong 
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E e (GeV) E e (GeV) 

Figure 9: Binwise distribution of the energy of e~~ in the final state e~e~$ for both signal-1 
and signal-2 (yfs = 500 GeV). The broken (— ■ ■—) brown lines correspond to the signal through 
production while the broken (—■■■—) green lines correspond to the signal events from e^-pair 
production. The background follows the notation ofFig\^ Each binsize is 10 GeV. (a) £ = 100/6 -1 , 
(b)£ = 500/r 1 . 

^-channel contribution from the photon exchange, as discussed earlier. This actually helps 
suppress the huge continuum background in the case of e~e~ ]fi final states, as shown in 
Table [2j The distribution in AR ee also shows a clear difference in behavior from the SM 
background, which can again be attributed to the fact that the angular distributions for the 
final state charged leptons in the signal events are markedly different from those of the SM 
events, and if one considers 

\&n\ = We- 

the distribution actually shows a peak for the signal at \An\ = 0, while for the SM background 
the signal peaks beyond \Arj\ > 1.2. 

Finally we present the energy distribution of the charged leptons for the two different 
final states in Fig. [9j The signal is again seen to stand out against the 5a fluctuations in 
the SM background for both signal-1 and signal-2 for the higher integrated luminosity of 
500 fb' 1 in Fig. MJo). The more interesting feature is seen for signal-1, where the energy 
distributions for the final two charged leptons will in principle be different. Assuming that 
it would be difficult to distinguish between the two leptons, we show the distribution by 
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Figure 10: Binwise distribution of missing energy for both signal-1 and signal-2 (y/s = 1 TeV). 
The background follows the notation of FigO Each binsize is 20 GeV. (a) £ = 100/6" 1 (hnal 
state e~e~E), (b) C = 500/&- 1 (final state e~e~E), (c) C = lO/b^ 1 (final state H~H~$). Here 
M A -- =500 GeV. For signal-1 m iL =250 GeV and for signal-2 mg L =500 GeV. 



taking the average of both the leptons energy distributions. However, in the case of signal-1 
the leptons coming from the 2-body decay of slepton will have an energy profile depending 
on the mass difference — M~o , while the lepton which comes from the 2-body decay of 

A will have an energy profile dominantly depending on ~ m i- For signal-2 however, 

both leptons have identical distributions which depend on the 3-body decay kinematics of 
the parent particle. 

In Figs. [TUJ dU and [12] we show the distributions for the total missing energy, invariant 
mass of the visible leptons and the energy profile of the final state leptons for the machine 
with -\/s = 1 TeV center-of-mass energy for the sample point B given in Table [H The mass 
for is taken as 500 GeV while the corresponding masses of the selectrons for signal-1 and 
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Figure 11: Binwise distribution of the invariant mass of the visible particles in the final state for 
both signal-1 and signal-2. Each binsize is 20 GeV (y/s = 1 TeV). The background follows the 
notation of Fig\^ (a) C = 100/6" 1 (final state e~e~$), (b) C = 500/6" 1 (final state e'e'E), (c) 
C = Wfb- 1 (final state Here y/s = 1 TeV and M^— =500 GeV. For signal-1 mg £ =250 

GeV and for signal-2 mg L =500 GeV. 



signal-2 are 250 GeV and 550 GeV respectively. The contribution to signal-1 from e^-pair 
production in this case is however very small (0.2 fb) and does have any significant effect 
on the distributions. In each of the plots we also show the distributions for the final states 
fi~ Here, in addition to distributions given for integrated luminosities £ = 100/6 -1 
and C = 500/6 -1 , we also include the case (c) £ = 10/6 x . The signal is quite large for 
heavier and the signal far overwhelms the background for the n~ or f,. It can be seen in 
all the Figs. [ToyTTl and [12] that we still have a large SM background and that an integrated 
luminosity of £ = 500/6 -1 is needed for a significant 5<r-signal for signal-2 for the e~e~f! 
final state. The SM background for the pTpTlfj final state is however still quite suppressed 
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Figure 12: Binwise distribution of the energy of e (u ). The background follows the notation of 
FigO Each binsize is 20 GeV (y/s = l TeV). (a) C = 100/6" 1 (final state eTe~$), (b) C = 500/&- 1 
(final state e^e'E), (c) C = 10/6" 1 (final state M">~# Here y/a = 1 TeV and M A -- =500 GeV. 
For signal-1 rrie L =250 GeV and for signal-2 m^ L =500 GeV. 



and one can observe the large rates for both signal-1 and signal-2 as compared to the 5a 
statistical fluctuations in the SM for £ = lOfb^ 1 . We preferred to show the distributions on 
a linear scale as opposed to the earlier plots, since the rates for both signal-1 and signal-2 
do not differ by a large factor. The distributions do not exhibit any additional new features 
as compared to the plots for the same kinematic variables at the \/s = 500 GeV machine. 
However a look at Table [2] suggests that a heavier A accessible at a yfs = 1 TeV machine 
will have a substantial SM background and the signal-background analysis merits discussion. 

3.2 The right-handed higgsino 

We briefly discuss the production and decay of the right chiral higgsino A^ - at the e~e~ 
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Figure 13: The production cross section as a function of the triplet higgsino mass (y/s = 1 TeV). 
The cross section for both A^~ and A^~. The solid (red) line corresponds to A^ _ production 
with the beam polarizations (+1,+1) while the broken (blue) line corresponds to A^~ production 
with the beam polarizations (-1,-1). The plots shown are for the sample point A given in TableUi 

collider, emphasizing distinguishing features with respect to A^~. In Fig. [13] we show the 
production cross section for its production and compare it with the production cross section 
of A^~. Note that the production of the right-chiral A^~ requires the electron beams to be 
dominantly right-polarized and can be quite large compared to the production cross section 
of the A^~, as shown in Fig. O The difference is mainly due to the difference in coupling 
of the neutralino with e — eL{R) which in this case reads, 



The production of the right-chiral A^~ will have a very clean signal as the SM background 
is completely reducible using polarization of the beams. We have already listed the SM 
background in Table [2] corresponding to different beam polarizations. The right-polarized 
beams effectively kill the SM background because the dominant contributions come from 
W-boson exchange which vanish when the electron beams are right-polarized. So the signal 
for the A^~ production is relatively background free. In fact, this suggests that the signals 
for the right-handed A^~ will be even more striking and hence much lower values of the 
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AL = 2 couplings can be probed. To highlight this fact and to illustrate the sensitivity of 
the signal to the AL coupling, we calculate the rates for different values of the AL coupling 
f ee . We plot the distributions in the missing energy, invariant mass of the visible leptons 
and the energy of the visible lepton in the final state for the e~e~]fi final state, coming from 
the e^-pair production and production, and the subsequent decays in Fig[TU for three 
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Figure 14: Binwise distribution of missing energy, invariant mass and energy of e~ for different 
values of the AL = 2 coupling f ee - The statistical fluctuations in the SM background is shown at 
3cr in black (dashed) lines. Each binsize is 10 GeV while the luminosity is taken as C = 500/6 -1 
(hnal state e~e~ Ifi). The green (—■■—) line stands for the independent contribution from e^-pair 
production for signal-1. Here y/s = 500 GeV and — =300 GeV. For signal- 1 shown in red 
lines, rrie R =150 GeV and for signal-2 shown in blue lines, =400 GeV. 



values of f ee =0.04, 0.05 and 0.1. We consider a luminosity of C = 500/6 -1 and show the 
3cr fluctuations in the SM background. All the figures clearly show that signal-1 stands out 
against the background for a AL coupling strength which can be as low as f ee = 0.05 and 
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even lower values can be accessed if we have final states. The green (— ■ — ) line 

stands for the independent contribution coming from e^-pair production for signal-1. For the 
polarization choice (+1,+1), the cross section becomes 15.8 fb. One can see that for f ee — > 

the distribution for signal-1 will converge on the green ( ) line as the contribution coming 

from the production becomes negligibly small. The plot shown for signal-2, however 
shows the direct reach for the coupling f ee as the selectron-pair production is kinematically 
disfavored for a right-selectron mass of 400 GeV. 

4 Summary and Conclusions 

The absence of AL = 2 processes in the SM provides an unique platform to search for 
beyond the SM scenarios. Such processes arise naturally at an e~e~ collider where they can 
be produced directly. LRSUSY models where R-parity is conserved, give rise to such doubly 
charged exotics which can be relatively light and accessible at future linear accelerators. In 
this work, we have studied in detail the production of a single doubly charged higgsino and its 
decay channels at the ILC operating in the e~e~ mode. We showed that the production cross 
sections for the A can be quite large if the beams are polarized and that enough statistics 
is available to study the signals arising from its subsequent decays against the dominant SM 
background. We have also shown how kinematic cuts can effectively limit this background 
and give a clear signal for LRSUSY at an e~e~ collider. We described how different final 
states (e~e~fi, fx~ 'fi~ f) compare against the dominant SM background. The use of polarized 
beams helps in producing doubly charged higgsinos of different chirality and hence proves 
to be an essential tool in distinguishing between the two chiral states of the A producing 
the same final states. We also find that the production of the right chiral higgsino A^ can 
be relatively background free and is more sensitive to lower values of the AL = 2 coupling. 
This is a major advantage over the production of these exotics at hadron colliders or in 
e + e~ collisions where ascertaining the chiral nature of A will be a non-trivial issue. A 
dominantly right (left)-polarized option for both the electron beams can achieve this quite 
easily at the e~e~ collider. As most studies of LRSUSY expect the doubly- charged higgsinos 
to be light, and even allow for the next to lightest supersymmetric particle (NLSP) to be 
the right-handed doubly charged higgsino [9], the resulting signals could be spectacular. 
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To summarize, in this work we have shown that an e~e~ collider can be an ideal machine 
to study the production of doubly charged higgsinos of the LRSUSY model. Through the 
choice of a few representative points in the parameter space of the model, we have shown 
signals coming from the production of A for two center-of-mass energies y/s = 500 GeV 
and y/s = 1 TeV. We have shown that through the final states £~£~f!we have clear signals 
for these exotics. We also find that when £ — /i, r the SM background is a six-body final state 
and hugely suppressed. We have also highlighted the advantage of using polarized beams in 
producing these exotics and suppressing the SM background. 
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5 Appendix 

In what follows, we use the conventions of [23]. 



5.1 Charginos 

In the supersymmetric sector of the model there are four doubly charged higgsinos A^~ , 5~£ + , A]/~ , 
and The mass terms relevant to these higgsinos are: 

C K = -M K ^A- L -5i + - M A __Ar^ + + h.c, (6) 

where in our notation = /i 3 . The model also has six singly-charged charginos, corre- 
sponding to \ L , X R , 4> u , 4> d , Af, and A%. 

The terms relevant to the masses of charginos in the Lagrangian are 

/ 




+ H.c. , (7) 
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where ip +T = (-iA£ , fc, $" u , A£ , A£) and ^ T = (-iA £ , -iA fl , 2d , 2u , £ £ , and 
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where we have taken, for simplification, /i^ = fi\. Here k u and /c<j are the bidoublet Higgs 
bosons vacuum expectation values (VEVs), v^ R and vs R are the triplet Higgs bosons VEVs, 
and M L ,M R the SU(2) L and SU(2) R gaugino masses, respectively. The chargino mass 
eigenstates x% are obtained by 



Xi = VijipJ, Xi = Uijtpj , i, j = 1, . . . 5, 
with V and U unitary matrices satisfying 

irxv 1 = m d . 



(9) 



(10) 



The diagonalizing matrices U* and V are obtained by computing the eigenvectors corre- 
sponding to the eigenvalues of XX^ and X^X, respectively. 



5.2 Neutralinos 



The model has eleven neutralinos, corresponding to Xz, Xz', Ab-l, 4>1 u , 4>2ui 0?d> 02d> 
K° R , 81, and 8 R . The terms relevant to the masses of neutralinos in the Lagrangian are 



C 



v = -\^ T Y^ + E.c. , 



(11) 
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where = ( _ iA o ^ _ tX R} _ zXb _ L} ^ ^ A o j 50 ) A°j, We omit 0°J T , which only 
mix with each other. The mixing matrix is: 
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The mass eigenstates are defined by 

X° = A^°(U = 1,2,...9), 
where iV is a unitary matrix chosen such that 

NZN T = Z D , 
and Z D is a diagonal matrix with non-negative entries. 



(12) 



(13) 



(14) 



5.3 Scalars 

The interactions between vector bosons and scalars arise from the kinetic energy term for the 
gauge bosons in the Lagrangian density. We denote by xl,r the Higgs fields before mixing, 
and by Ul,r the Higgs fields after mixing. The Higgs scalar fields are defined as: 
Doubly Charged Fields 



x 



++T 

R — 



Vr 

X L 
±±T 

Vl 



(a£ + s- r -*) , x- R T = (A++* 6JT) 

(Ar si-*) , x- L - T = (a^ 5- L -) 



(15) 
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Singly Charged Fields 



* +T = (a£ 5 L * 0- 0+ 2 ; <t>t d a+ *-) , 

= (Ar <fe tfj 2u 0i* A+* ^) , 



1/ 



(ff* fl* fl* Ht Hf Hf Gf Gf) (16) 



Neutral Fields 



= (h Al H Sl Hl d H^ u H\ u Hl d H Ar H 5l ^j , 
yf ee (flj flj tf 3 ° iJ 4 ° Hi Hi H° Hi) , 

™0T _ / ^0 J3 ^0 „ \ 

x p — \ z &l z $l z ld z 2u z lu z 2d Z A R Z$ R j , 

y f ee (A? A° A° A° A° A° G? G% (17) 

where the indeces "s" and "p" stand for scalar and pseudoscalar, respectively. There are two 
charged Goldstone bosons for the left-handed and the right-handed charged vector bosons, 
and two neutral Goldstone bosons for the Z L and Z R bosons. They have zero mass. We 
define them to be the 7-th and 8-th component of H ± and A , in order to simplify the 
summation convention. The mass matrices M are real and symmetric and diagonalized by 
orthogonal matrices R defines as: 

= «,..., 

K% = df>« («•?!. 0,0) (18) 

where m sl is the mass of the lightest Higgs scalar. Diagonalizing the scalar mass matrices, 
we introduce new fields by: 

«°„ = * = K),,4, do) 
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5.4 Couplings 

Finally, we list the all couplings involving doubly charged higgsinos relevant for the produc- 
tion and decay: 



ki + gN k3 jP R (20) 

- gN k ^ P L (21) 

(22) 
(23) 
(24) 
(25) 

where is the matrix which diagonalizes the mass matrix for the singly-charged Higgs 
bosons. C is the charge conjugation operator while Pl and Pr are the chirality projection 
operators (1 =F7 5 )/2. 
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